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Abstract 

Data on KeA, decays allow one to extract experimental information on the elastic 
TTTT scattering amplitude near threshold, and to confront the outcome of the analysis 
with predictions made in the framework of QCD. These predictions concern an isospin 
symmetric world, while experiments are carried out in the real world, where isospin 
breaking effects - generated by electromagnetic interactions and by the mass difference 
of the up and down quarks - are always present. We discuss the corrections required 
to account for these, so that a meaningful comparison with the predictions becomes 
possible. In particular, we note that there is a spectacular isospin breaking effect in 
decays. Once it is taken into account, the previous discrepancy between NA48/2 
data on Ke4, decays and the prediction of vrvr scattering lengths disappears. 
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1 Introduction 



Chiral perturbation theory (ChPT) [1], combined with Roy equations, allows one 
to make very precise predictions for the values of the threshold parameters in 
elastic vrvr scattering [2] - for a status report, see e.g. the contribution of one of 
us at KAON'07 [3]. Several experiments have measured the tttt interaction at low 
energy with a precision such that it is now possible to confront these predictions 
with data: i) Ti'^Ti^e^Vf. decays [4-8], ii) the pionium lifetime, measured by 

the DIRAC collaboration [9], iii) the cusp effect in K 3n decays, investigated 
by the NA48/2 and KTeV collaborations [8,10]. The experiments performed by 
the NA48/2 collaboration have generated an impressive data basis, as a result of 
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which the matrix elements of K — s> Svr and of K^a decays can be determined with 
an unprecedented accuracy. 

The theoretical predictions and the measurements are performed in two dif- 
ferent settings: the predictions concern pure QCD, in the isospin symmetry limit 
mu = rnd, with photons absent - a paradise world. To be more precise, the con- 
vention is to choose the quark masses and the renormalization group invariant 
scale of QCD such that the pion and the kaon masses coincide with the values 
of the charged ones, and the pion decay constant is = 92.4 MeV. [We do 
not specify the masses of the heavy quarks, because in the present context, their 
precise values do not matter.] On the other hand, experiments are all carried out 
in the presence of isospin breaking effects, generated by real and virtual photons, 
and by the mass difference of the up and down quarks: this is the real world, de- 
scribed by the Standard Model. We are thus faced with the problem to find the 
relation between quantities measured in the real world, where isospin breaking 
effects are always present, and the predictions made in the paradise world. It is 
the aim of the present article to provide this relation. 

In early experiments [4], the effects of real and virtual photons were estimated 
by considering a simplified model for the weak interactions, and taking into ac- 
count photon effects through minimal coupling, working at lowest nontrivial order 
in ctQED- In analyses of NA48/2 data before spring 2007, real and virtual photon 
effects were treated in a factorized manner, applying the Coulomb factor and 
using the program PHOTOS [11] - see Ref. [5] for details. It is clear that both 
approaches missed the effects generated by the pion and kaon mass differences, 
and by the quark mass difference rrid — rriu. It turned out that these effects are 
quite spectacular [12-15]: when taken into account, all previous discrepancies 
between data and prediction for the scattering lengths disappear [3,6-8jl|. 

The outline of the paper is as follows. In section [2|, we recall how K^a data are 
analysed in the isospin symmetry limit, based on Watson's theorem, while section 
|3] describes the framework in which we take isospin breaking effects into account. 
In section H] we investigate - for the case of (Lorentz) scalar pion form factors 
- the effect of isospin breaking mass differences using unitarity and analyticity 
alone. We illustrate the outcome of this investigation in section[5]with two explicit 
examples in the framework of Quantum Field Theory, before we return to K^a 
decays in section [6l where we display the result for the changes in the phase of the 
form factors, and for the phase-removed form factors themselves. These results 
allow us in section [7] to perform fits to Ke4, data, based on numerical solutions to 
Roy equations, and to determine in this manner experimental values for the tttt 
scattering lengths. A comparison with related work available in the literature is 
provided in the following section [8], while a summary and concluding remarks are 

"'^In fact, the relevant expressions for the pertinent corrections are contained ah'eady in the 
early works of Cuplov and Nehme [16], but they went, unfortunately, largely unnoticed in the 
literature. 
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given in section [91 Appendix [B] contains proofs of the general statements made 
on the low-energy structure of phases and form factors in the main text. Finally, 
Appendix [C] contains material needed for the calculations of form factors in a 
non-relativistic framework. 



2 Ke4 decays: isospin symmetry limit 

To set notation and to explain the manner in which Kg^ decays allow one to mea- 
sure TTTT phase shifts, let us consider here the decays in the isospin symmetry limit 
iTT'u = ''^d^ttQED =0. The matrix element for K~^{p) '^^{Pi)'^~{P2)^^{Pe)^e{Pu) 
is 

T = ^v:MP.)ni - iMpeW, - A,), (1) 

where the last factor denotes hadronic matrix elements of the strangeness chang- 
ing (vector and axial vector) currents, 

V^- Af, = (7r+(pi)7r~(p2) out | {s-f^u - s'jf.-f^u) \ K^{p)) ■ (2) 

In the following, we concentrate on the matrix element of the axial vector current, 
because it carries information on the vrvr final state interactions and, in particular, 
on the TTTT phase shifts. One decomposes into Lorentz scalars, 

= -^TT [(^1 + ^^2)^^ + (Pi - P2)^G + (pe + • (3) 

The form factors F, G, R are holomorphic functions of the three variables 

S = {pi+ P2f ,t = {pi-pf ,U = {p2 - pf . (4) 

Sometimes, it is useful to use instead 

S= {Pl+ P2Y , Si = {pe + Puf , COS 6^^ , (5) 

where 9-,^ is the angle of the vr"*" in the CM system of the two charged pions, 
with respect to the dipion line of flight in the rest system of the kaon [17]. In the 
isospin symmetry limit, one identifies the vrvr phase shifts in the matrix element in 
a standard manner, by performing a partial wave expansion, and using unitarity 
and analyticity, although, in the present case, this is a slightly intricate endeavour 
[18]. It is useful to introduce a particular combination of form factors (we omit 
isospin indices). 
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Here, a = a/1 — 4M2 / s, and X{x,y,z) is the triangle function. The form factor 
Fi has a simple partial wave expansion, 

F, = Y,Pkicose^)h{s,s,). (7) 

fe>0 

For fixed si, the amplitudes fk are holomorphicH in Cl/j(4M^). In the elastic 
region, the form factors fk carry the tttt phase shifts in the pertinent isospin 
channel [18]. For /q and /i, the relation is 

U = e2^^°/o-,/+ = e2*^Vr; [4M2,16M^2], (8) 

where 5o (^i) denotes the phase shift of the isospin zero S-wave (isospin one 
P-wave), and fnifn) stands for the form factor evaluated above (below) the cut, 

= fn{s±te,s,). (9) 

The phase-removed form factor can be Taylor expanded at threshold. For the 
quantity /o, the expansion reads 

e-''^fo'~ = co + c,q' + 0{q'), 9^ = ^-1, (10) 

with coefficients Cj that depend on Sf. Because the modulus squared enters 
the decay rate, one can measure the phase shift difference 6q — 6i in K^i decay 
experiments. In the remaining part of this article, we investigate the manner in 
which the relations Eqs. (IHfTOll are modified in the presence of isospin symmetry 
breaking effects, and how the analysis of Ke4 data must be modified in order to 
determine the tttt phase shifts. 

Before proceeding, let us consider the form factors at one loop in chiral per- 
turbation theory and verify that Fi indeed has the behaviour just discussed. For 
this, we recall the pertinent effective Lagrangian 

C2 = ^{D^UDf'U^ + 2BoMiU + U^)), (11) 

where the covariant derivative D^U contains the external vector and axial vector 
currents, and A4 = diag(m, m, m^). Some of the graphs that contribute at tree- 
level and at one loop are displayed in figure [H The full result is [19] 

Mk 



with 



/o(., s,) = -J^ {1 + A{s) + H{s, s,) + 0{p')} , (12) 



A{s) = ±-^{2s-Ml)J{s), 

167r2J(s) = a ( In ^— ^ + ZTT ) + 2 , s>AMI. (13) 
V 1 + o- / 



^See Appendix \K\ for our notation of the various sets in the complex plane. 
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a) b) c) 

Figure 1: Some of the graphs that contribute to the matrix element of the axial 
current at tree and one-loop order. The filled vertex indicates that the axial 
current also couples to a single kaon line. That graph contributes to the form 
factor R. There are many additional graphs at one-loop order, not displayed in 
the figure. 

Here, (-^o) denotes the pion mass (pion decay constant), at leading order in 
the chiral expansion. The quantity H{s, se) is real in the interval of elastic tttt 
scattering. It is now seen that /o indeed has the property ([8]) at this order in the 
low-energy expansion, with 

This is the phase shift of the isospin zero S-wave, in tree approximation. The 
amplitude /i has a very similar structure [19], containing the phase shift of the 
isospin one P-wave, again in tree approximation. Remark: To get the result 
Eq. (|T4l) . we have used an expansion around the chiral limit m„ = = = 0, 
as a result of which the pertinent pion decay constant Fq in this limit appears. 
To the order considered here, we may replace Fq by the pion decay constant F in 
the chiral limit = = 0, ^ 0. It seems to us that this is a more natural 
choice when discussing the vrvr phase shift, and we will, therefore, use Fq F in 
numerical analyses in the rest of this article. 



3 Isospin breaking: the framework 

The relation ([8]) holds in the isospin symmetry limit. On the other hand, K^^ 
decays happen to occur in the real world where 7^ niu, oqed 7^ 0, as a result 
of which the relations Eq. ([8]) do not hold anymore. Here, we discuss how this 
situation can be modelled in view of the already published analyses of Kf,^ decays. 

In early experiments [4], the effects of real and virtual photons were esti- 
mated by considering a simplified model for the weak interactions, and taking 
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into account photon effects through minimal couphng, working at lowest non- 
trivial order in oqed- In analyses of NA48/2 data before summer 2007, real and 
virtual photon effects were treated in a factorized manner, applying the Coulomb 
factor and using the program PHOTOS [11] - see Ref. [5] for details. It is clear 
that both approaches missed the effects generated by the pion and kaon mass 
differences, and by the quark mass difference m^ — m^j. These must thus be taken 
into account separately. We assume that they can be evaluated in factorized form 
as well, and write symbolically in case of the NA48/2 analysis 

Full isospin breaking effects = Coulomb factor x PHOTOS x mass effects 

As a practical way of proceeding, which should catch the main effects, we propose 
to correct also earlier analyses [4] with the last factor. 

We now discuss the manner in which mass effects may be evaluated. As 
the correction turns out to be small (although not negligible), a perturbative 
method is appropriate. In the following, we use effective field theory techniques 
to perform the calculation. As real and virtual photons have already been taken 
into account, we need a framework that accounts for mass effects only. An obvious 
candidate is the chiral lagrangian itself, in the absence of real photons. This can 
be achieved by modifying the lagrangian £2: one adapts the quark mass matrix, 
Ai — > diag(mu, m^, mg), and adds mass breaking terms of electromagnetic origin 



where C is a low-energy constant that breaks the isospin symmetry of the meson 
masses, such that Mj^ 7^ Mt^o,Mk 7^ M^o. The symbol e stands for the electric 
charge. At order p^, there will be additional terms, classified in Refs. [20,21]. 
The effect of the replacement ( !T5|) is twofold: first, because the pion masses split, 
the loop contributions generated by the diagrams displayed in Fig. [T]d),c) have a 
different threshold. Second, in addition to the graphs displayed in figure [H there 
is a new contribution shown in figure [2J the kaon interacts with the axial current 
to generate a Ti^rj intermediate state. Because 7^ m^, the 77 can transform 
back into a neutral pion, that then re-scatters with the second neutral pion into 
a charged pion pair. We perform later in this article a quantitative analysis of 
these effects. 

In summary, we propose the following framework to purify published phase 
shifts from mass effects. 

i) We assume that the published phase shifts correspond to the ones obtained 
in a world defined by the lagrangian f|T5l) . 



[20], 



C2^C2 + C{QUQU^) , Q = ^diag(2, -1, -1) , 



(15) 



e''^" f^,k = 0,1- s>4Ml 



(16) 



We refer in the following to ipk as "measured phase shifts" . 
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Figure 2: One-loop graph in the Kei decay amphtude that emerges due to the 
7r°?7 mixing. There is no counterpart of this graph in the scalar form factor. 

ii) In order to get into contact with the isospin symmetric phase shifts 6^, we 
note that 

h = ^k- ii'k -Sk),k = 0,1. (17) 

The differences ipk — 6^ on the right-hand side are small, and can be calcu- 
lated in the effective field theory framework outlined above. After subtract- 
ing these from the measured phase shifts ipk, one gets 6k, which can then be 
confronted with predictions in the framework of ChPT. 

While this procedure does not present a complete and full analysis of radiative 
corrections in K^.^ decays, it allows one to purify published phase shifts from mass 
effects and thus to hopefully retain the main effects of isospin breaking. [One may 
envisage a more ambitious procedure [22], by working out the relevant matrix 
elements in the framework of ChPT including photons and leptons [23] , and then 
constructing a new event generator, to be used in a new analysis of K^i decay^. 
Eventually, such an analysis might lead to an improved algorithm. However, we 
consider this to be a long term project.] 

The framework proposed here allows one to investigate the strictures imposed 
on the form factors by unitarity and analyticity, even without relying on any spe- 
cific details of the lagrangian Eq. f|T5l) . We find it instructive to discuss this fact 
in some detail in the following section - before performing the explicit evaluation 
of isospin breaking effects in decays - because it illustrates that, qualitatively, 
the effects we are finding at the end in those decays do not rely on a specific la- 
grangian framework, but are present in any theory that breaks isospin symmetry 
and incorporates unitarity and analyticity. [On the other hand, to work out the 
effects in a quantitative manner, a specific underlying theory is needed.] Readers 
who are not interested in this general setting may wish to skip the following two 

This would also take care of the Coulomb phase, whose effect is suggested to be substantial 
in Refs. [13, 14, 16]. While the Coulomb phase acts in the same direction as the mass effects 
considered here ~ it increases the difference ipo — ipi — {5o — di), see section[7]- it is not clear which 
part thereof is already included in PHOTOS. We, therefore, prefer to stick to the procedure 
proposed here, because these arc corrections that arc definitely not included in PHOTOS. 
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sections and continue directly with section [6l where we take up the case of Kg^ 
decays again, and where we calculate the corrections ipo — 60 and -i/^i — Si, and 
investigate their effect on the tttt scattering lengths. 



4 Watson's theorem and all that 

Here, we discuss Watson's theorem and its modification in the isospin breaking 
case by using analyticity and unitarity arguments alone. As already mentioned, 
it is not quite straightforward to work out unitarity and analyticity constraints 
on the axial current matrix element Eq. ([3]) which is relevant here, because this 
matrix element describes the scattering process K + axial current — > vr + vr, which 
has a rather complicated structure. On the other hand, the physical effects 
generated by isospin breaking interactions are also present in simpler matrix 
elements like the (Lorentz) scalar form factor of the pion. We, therefore, illustrate 
the basic facts in this simpler setting. It will become obvious that the same line 
of reasoning could be carried over in the same general framework to the Kg^ form 
factors, although with much more labor. 

To set the framework, we consider matrix elements of a hermitian Lorentz 
scalar current which is taken to be isoscalar in the absence of isospin break- 
ing interactions. To simplify the analysis further, we consider the case where the 
isospin symmetric theory corresponds to the standard SU (2) SU (2) ^ — > SU{2) 
scenario of spontaneously broken two-flavour QCD, with pions only. Emission of 
real photons is excluded, as a result of which all matrix elements are infrared 
finite. 

As isospin symmetry is assumed to be broken, it is convenient to use state 
vectors that are labelled by the physical pion states. The two form factors are 



In the isospin symmetry limit, one has = Fq in the Condon-Shortley phase 
convention used here. These form factors are boundary values of functions Fk{z) 
which are assumed to be i) holomorphic in Cji{4:M^o), and ii) real on the real 
axis for s < 4M^o- For further specifications of Fj. needed to arrive at the results 
described below, see Appendix |Bl 

We also need the elastic tttt — > vrvr scattering matrix elements and consider the 
following three channels in vt^tt^ iiH'^ : {ab; cd) = (00; 00), (H — ; 00), (H — ; H — ). 
We denote the scattering matrix elements by Ti{s,t), where i = 0,x,c are the 
channel labels in the above order, and s, t stand for the Mandelstam variables. 
The partial-wave expansion reads 



(0|j(0)|vr+(pi)7r-(p2);in) 






00 



Ti{s,t) = 32n /ii(s) + ^Pi(cos^9)/i[(s) 



(19) 



1=1 
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Here, 'd stands for the scattering angle, and hi{s) denote the pertinent S-waves, 
which are needed in the following. They are boundary values of functions hi{z) 
which are assumed to be i) holomorphic in Clr{4:M^o), and ii) real on the real 
axis for < s < 4M^o. For further specifications of hk needed to arrive at the 
results described below, see Appendix El 

The unitarity conditions for the form factors and for the partial waves read 
on the upper rim of the cut 

ImF(s) = T(s)p(s)F*(s) , 

ImT(s) = T{s)p{s)T*{s) ; AM^o < s < IGM^o . (20) 
We have used the matrix notation 

F 







he 










ho) 



^^'^"^ ao(.)^(.-4M2,) ) ' ^^^> 



together with 



, .„M^,1-^^. (22) 



The quantity a{z) is holomorphic in Cl/j(4M^). On the upper rim of the right- 
hand cut we take a{z) to be positive and real, and the analytic continuation 
thereof elsewhere in the complex z-plane. Analogous statements hold for ao{z). 

4.1 Isospin symmetry limit 

In the isospin symmetry limit, one has 

{K ,h,,ho) = l (^2 + 2to , 2t2 - 2to , 4t2 + 2to) • (23) 



The quantity tj denotes the partial wave amplitude with angular momentum zero 
and isospin /. With Fc = Fq, the unitarity relations Eq. fl20|) decouple. 



ImF, = toaF:, 

Imtj = 1 = 0,2. (24) 

We analyse in Appendix [B] the singularity structure of the partial waves and of 
the form factor which follow from these relations [24]. The result is as follows: 
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i) The form factor develops a square root singularity at the threshold z = 4:M^, 

Fc{z) = A{z) +iaB{z) , (25) 

where A, B are meromorphic in Clr(16M|), and real on the real axis in the 
interval < s < 16M^. 

ii) Define the phase 

e{z) = arctan^^^. (26) 
^ ' A{z) ^ ' 

In the elastic interval, 9 coincides with the phase shift 5q of the partial 
wave to ^ this is Watson's theorem. The above definition of 9 holds also 
for complex values of 2;. As a result, 9 is holomorphic in the shaded region 
shown in Fig. [HI cut along the positive real axis for s > 4M^. [Although 9{z) 
is complex in general, we keep calling this quantity phase for simplicity.] 

iii) In analogy to Eq. ( ITOl) . we define the phase-removed form factor 

F^{z) = e-'^(^)F,(z) = A (^1 + ^ j . (27) 

Fc is holomorphic in the shaded region shown in figure [3l and real on the 
real axis in that region. It can therefore be Taylor expanded at threshold, 

F. = 5^e,(;.-4M,2)% (28) 

with real coefficients Cj, cf. with Eq. (fTOl) . 

Analogous facts were (implicitly and explicitly) used for the axial form factors fk 
in Eq. (jTj) in all previous analyses of K^.^ decays. 

4m! 



Figure 3: Region of analyticity in the complex 2;-plane, isospin symmetric case. 
The phase-removed form factor Fc{z) is holomorphic in the shaded region. The 
symbol Zimi. denotes the inelastic threshold z = 16M^. 
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Figure 4: Region of analyticity, isospin broken case. The phase-removed form 
factors Fk are holomorphic in the (cut) shaded region. The symbol Zinei. denotes 
the inelastic threshold z = 16M^o- 



4.2 Broken isospin symmetry 

We now consider the case where the underlying theory contains isospin breaking 
interactions. To render the discussion simple, we assume that these terms are 
reasonably small, such that the isospin symmetric values of the form factors and 
of the scattering matrix are only slightly modified. This can always be arranged 
by adjusting the relevant parameters in the lagrangian, and allows us to discuss 
the changes induced in a simpler manner. We analyse the singularity structure 
of the form factors Fc^o in Eq. f|T8l) in the general case in Appendix [B] as well, 
and we find that the four statements i) to iv) made above are modified in the 
following manner: 

i') The form factors have the threshold behaviour 

Fk = Ak + iaBk + ia^Ck + crcroDk ; /c = c, , (29) 

with coefficients Ak, . . . ,Dk that are meromorphic in CLij(16M^o), and real 
on the real axis in the interval < s < 16M^o- 

ii') Let 

= arctan ^'fff ' 7°'^^^';' ■ (30) 

Ak{z) + a{z)ao{z)Dk{z) 

and introduce the phase-removed form factors Fk, 

Fk{z) = e-'''^^'^Fk{z). (31) 

The quantities Fk are holomorphic in the (cut) shaded region shown in Fig. 
m and real on the real axis in the intervals < s < 4M^o and 4M^ < s < 
1QM\. 

iii') In contrast to the isospin symmetric case one finds that 

a) the vrvr scattering amplitude does not fully determine the phases 6k 
below the inelastic region, 

b) 9c does not vanish at the threshold s = 4M^ , 

e^^O^s^AM^. (32) 
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Figure 5: Bubble diagrams, contributing to the neutral scalar form factor in the 
non-relativistic effective theory. Solid and dashed lines correspond to charged 
and neutral pions, respectively. 

iv') Again in contrast to the isospin symmetric case, the form factor F(. develops 
a square root singularity ai z = 4M^, 

F,{s) = eo + eiq + 0{q^) ,s\AMl; (33) 

with real coefficients . 

This result shows that the phase-removed form factor cannot be expanded 
in a Taylor series at the threshold s = 4M^. [One may be tempted to get rid 
of the problem by defining a form factor where the pertinent Omnes factor is 
removed. In the present case, this is not what one wants to do, because the aim 
is to measure the phase shift, which is needed to evaluate the Omnes factor.] 

5 Explicit examples 

We find it useful to provide in this section examples of form factors that illustrate 
the analytic properties worked out above. 

5.1 Scalar form factor in a non-relativistic effective theory 

In order to describe the behaviour of the scalar form factor in the low-energy 
region = s/{4:M^) - 1 < 1, = - M^o < we use the framework 
of non-relativistic effective field theory (NREFT) as developed in Refs. [25-28]. 
This formulation is especially convenient to study the singularity structure near 
threshold, because the locations of the low-energy singularities coincide with 
that in the relativistic QFT to all orders in the low-energy expansion. A similar 
non-relativistic framework has recently been used to study K ^ Sir decays in 
the vicinity of cusps [25-27]. The effects that we are addressing in this section 
have the same physical origin as the cusps in K ^ Sir decays. They emerge, 
because the final state interactions involve both, tt^tt" and tt'^tt^ pairs, with 
different masses M,r 7^ M^o. We display the pertinent non-relativistic lagrangian 
in Appendix [Cl 

The scalar form factors Fc{s) and Fq{s), defined in Eq. (lC.16p . are given 
by a sum of bubble diagrams, see Fig. O The infinite series of these bubble 
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diagrams can be explicitly re-summed. The resulting expressions for Fc{s), Fo{s) 
are boundary values of meromorphic functions Fc{z), Fq{z) at 2 — s + iO, with 



1 — 2io-odo — Aiadc — ^aaox 



Fo{z) = — —, (34) 

1 — 2iaodo — Aiadc — 2aaox 

where di, i = c,x,0 and fk, k = c,0 denote polynomials of first order in z, and 
X = 4:{dcdo — dl). The coefficients of these polynomials are expressed through 
various non-relativistic couplings which, in turn, are determined by performing 
the matching to the underlying relativistic theory (see, e.g. [25-28]). In particular, 
the polynomials di contain 4-pion non-relativistic couplings and are expressed in 
terms of the effective-range expansion parameters for the tttt scattering through 
the matching of the tttt amplitudes at threshold. Similarly, the polynomials fk 
are determined from the matching to the relativistic form factors, expanded at 
threshold. 

Further, the quantities di, fk contain isospin-breaking corrections. We write 
di = di + d'i (and similarly for fk), where di stands for di calculated in the isospin 
limit. In analogy with Eq. (1231) . di can be expressed through two first-order 
polynomials ^0,^^25 corresponding to total isospin J = 0, 2 , 

{dc, 4, do) = ^ (^2 + 2^0, 2^2 - 2^0, 4^2 + "^Vq) . (35) 

Here, 

z - 4M2 

viiz) = ai+ ^j^^"" ain, 1 = 0,2, (36) 

and aijaiTi stand for the scattering length/effective range parameter in the 
isospin limit, see Ref. [25]. 

Providing similar explicit expressions for the isospin-breaking corrections in di 
and fk is not possible in general. The form of these corrections is not universal and 
the calculations should be performed within a particular underlying relativistic 
theory. An example of calculations in ChPT at one loop is considered in the 
following subsection. 

We add several comments concerning the structure of the two form factors. 

i) As already mentioned, the form factors Fj^z), F^i^z) are meromorphic func- 
tions in the cut plane displayed in Fig. [6^. The non-relativistic domain in 
the complex plane is defined as a strip surrounding the physical cut - it is 
indicated with the shaded region in the figure. It extends slightly below and 
above of the physical branch point, well below the first inelastic threshold. 
The maximal distance to the branch point in the non-relativistic region is 
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Figure 6: The analyticity domain of the form factor -Fc(^) (left panel) and the 
phase-removed form factor Fc{z) (right panel). The shaded area denotes the 
non-relativistic region and the filled circles indicate the positions of the poles on 
the first Riemann sheet. Both the left-hand cut and the poles lie outside the 
non-relativistic domain. 



set by the mass scale M^^. NREFT makes sense only in the non-relativistic 
region - the cut along the negative 2;-axis, as well as the distant poles on the 
first sheet should be regarded as artefacts of the non-relativistic treatment. 

ii) The form factors have the representation (!29|) . Constructing the vrvr amph- 
tudes in the same manner as the form factors, one finds that the unitarity 
relation given in Eq. fl20l) indeed is satisfied. 

iii) In order to obtain the phase-removed form factor, we define the phase as in 

Eq. dsni), 



arctan 



P-Q 
l + PQ 

2aodo + Aadc 



1 - 2aaox 
The phase-removed form factor 



Q = 2ao(^do+^jd^ . (37) 



2crcroX 



p2 



1/2 



(3^ 



is holomorphic in the (cut) shaded region shown in Fig. [6b. This differs from 
the isospin symmetric case, where the cut from 4Af^o to 4M^ is absent. 



15 



Figure 7: The two-loop diagram producing the hnear term ~ a in the phase- 
removed form factor, see Eq. 



In conformity with Eq. fl33p . the expansion of Fc{z) at threshold contains 
even and odd powers of the variable q, 

F,{z) = F,(4M2)(l+eig + e2g2 + 0(g=')), s\4M2 



ei 



-84'M-iAy^ + 0(A^/2), (39) 



where is evaluated at the threshold z = 4M^. Using Eqs. (1551) and (15^ . 
we get 

ei = (ao - a^)^ + 0(A^/2) ^ -0.004 . (40) 

At lowest order, this term is generated by the two-loop diagram displayed 
in Fig. [71 Numerically, the effect is tiny. 

iv) Expanding the phase in powers of a, ao gives 

e,{z) = Aad, - 2ao4y + 0({a, ^o}') • (41) 

Jc 

This expansion is very useful in the context of ChPT, since the generic 
expansion parameter here is {cr, (Tq} x c^c.x.O) with (ic,a;,o = OijP'). Thus, the 
neglected terms in Eq. (14T!) are of order and higher. Furthermore, in 
order to calculate the phase at 0{p'^) (two loops), the ratio /o//c should be 
evaluated at O(p^) and the matching of the polynomials d^d^ should be 
performed at 0{p'^) (at this order, di{z) are polynomials of second order in 
z). To this end, it suffices to use the one-loop result both for the form factor 
and for the tttt scattering amplitudes. 

v) As mentioned earlier, the phase 9c is not determined by the tttt scattering 
amplitude alone. Indeed, 6c depends on the ratio fo/ fc, which can take 
any value in the absence of isospin symmetry. In addition, this ratio con- 
tains low-energy constants (LECs) which do not occur in the tttt scattering 
amplitude, as will be illustrated in subsection 15.21 

vi) The standard way to define the phase of the form factor on the real axis in 
the interval s > 4M is 

lmFc(s + iO) , , 

09c = arctan — — — . (42) 

^ ReFc(s + iO) ^ ' 
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Figure 8: The phases (fc and as well as their difference. The cusp in the phase 
^Pc at the charged threshold s = 4M^ is clearly visible. The isospin-symmetric 
parts di have been determined from the matching condition, Eqs. (^5^ and 
with the scattering lengths and effective radii taken from Ref. [2]. The isospin 
breaking corrections d'^ were evaluated in ChPT at O(p^), see Eq. fl50l) . The ratio 
/o//c was set to one. 



The quantity ipc does not vanish at the charged threshold s = 4M^ and has 
a cusp there. This is illustrated in Fig. [8], where the phase ipc is plotted. For 
comparison, the isospin-symmetric phase 6o and the difference (fc — are 
displayed as well. The phases 6c and (pc are identical on the real axis above 
the charged threshold, 

Ms) = 9c{s) , s > ml . (43) 

On the other hand, they differ e.g. in the interval s G [4M^o,4M^], because 
6c{s) becomes complex there. 

5.2 Form factor of pion: chiral expansion 

We calculate the scalar form factor of the pion in the SU{2) version of ChPT. 
The Lagrangian is given by [20] 

£ = £2 + A , (44) 

where 

A = —{d^Udf^U^ + 2Bs{U + U^)) + CiQUQU^) . (45) 
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In the above expression we use standard notations. In particular, F denotes the 
pion decay constant in the chiral hmit, U is the pion field matrix, the quantity 
B is related to the quark condensate, s = Ai + s' is the scalar source, Ai and 
Q stand for the quark mass and charge matrices in the SU{2) case, respectively, 
and the coupling C is proportional to A^^ at lowest order in the chiral expansion. 
In accordance to the discussion in section [3l we do not include virtual photons. 

We do not display the order Lagrangian £4 explicitly. It is given, e.g., 
in Refs. [20]. Note that, since we do not include virtual photons, the ultraviolet 
divergences of the "electromagnetic" LECs ki from Refs. [20] change in an obvious 
manner. 

In the calculations, the scalar current j{x) = {2B)~^ {u{x)u{x) + d{x)d{x)) 
was used. The scalar form factors of the pion at one loop can be calculated in a 
standard manner and are given by the following expressions. 



Fc{s) = l + ^^is + AA^)J{s) + is-M^o)Ms) 



327r2F2 



{s + 4 A^)(L + 1) - (s - M^o){Lo + 1) + M^oLc 



Fo{s) = l + ^^i^2{s-Ml,)J{s) + MlJo{s] 

+ {-2(s - Ml,){L + 1) - Ml,{L, + 1) + 2MlL - Ml^L, 

+ -jfll + j-J. + e'K + Oip'). (46) 



Here 



L = ln^, Lo = ln^, (47) 

and the loop function J(s) is displayed in Eq. f|T3|) [Jq is obtained from it by the 
replacement cr — (Tq]. Further, 

20 4 
K = --{kl + k;) + -{5kl + 23kl + ^) + 8kl, 

20 4 
^5 = --ikl + k;)-2i-2k; + kl) + -i5kl + 5kl + kJ;). (48) 

In these expressions and /c[ denote scale-dependent renormalized LECs, the k^ 
adapted to the framework considered here (no photon loops). The quantity /i is 
the scale of dimensional regularization in ChPT. 
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The scalar form factor of the pion in the presence of isospin breaking have 
been evaluated in Ref . [29] . Our expressions agree with those of Ref . [29] up to 
contributions of virtual photons and terms of order e^, e^(mrf — m„) and {rrid — 

The phase of the charged form factor at one loop, extracted from Eq. (H6l) . is 
given by 

^^(^) = {(4^- + + ^^'O^o} . (49) 

It is seen that 6c indeed does not vanish at the threshold s = 4M^. 

Next we determine the non-relativistic couplings di, fk, performing the match- 
ing to ChPT at O(p^). We start from d^, which can be directly read off from the 
tree-level S-wave vrvr scattering amplitudes: 

^o(^) = (50) 

The quantities fk at 0{p'^) can be determined by performing the matching to the 
one-loop expressions for the form factors, given in Eq. (146|) . For the ratio /o//c 
one gets 



^ = H^ + H^{s-AMl)+0{{s-AMlf) 

Jc 



Ho = l+^|^ + e'(C;-Cn+0(A^,p'), 

TV 

We note that in the case of isospin violation the ratio /o//c is not determined 
only by the LECs that enter the vrvr amplitude. For example, even in the chiral 
limit the ratio /o//c contains the LEG kg, which is absent from the vrvr amplitudes 
in this limit. Consequently, in the case of broken isospin, the phase of the form 
factors is not determined by the vrvr scattering amplitudes alone. 

It is easy to check that the expression for the phase Eq. P9l) can be directly 
obtained from Eq. fHT]) by using the result of the tree-level matching, given in 
Eq. fl50|) . At this order, one may take /o//c = 1 + O(p^). 

In order to check the convergence of the chiral expansion, we have evaluated 
the phase of the charged form factor at two loops. We do not display the final 
result here. As expected, the next-to-leading correction lies within approximately 
30% of the leading-order result, see section [7] for a more quantitative result. 
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6 decays: isospin breaking 



We now come back to the analysis of Ke4 decays. We simply need to repeat the 
above analysis, this time for the matrix element Eq. ([3]) of the axial current. 

Working out the contributions from diagrams Figs. [1] and [21 one finds that 
the phase ipQ of the form factor /q [see Eq. (fT6|) ] becomes in the elastic region 
4M2 < s < 16M% 

= 32^ {^^^^ + + ~ + Ir) ^4 + ^^^'^ ' ^^^^ 

with 

rris — m , , 

R = — , 53 

nid - ruu 

and F the pion decay constant in the SU(2) chiral limit. The phase of the form 
factor /i does not contain at this order any isospin breaking effects [these are 
generated by tt^tt^ intermediate states, which cannot couple to P- waves], as a 
result of which one has 



i^i = 5i (54) 

at this order in the low-energy expansion. The one-loop expressions for the form 
factors F, G given in Refs. [16] contain the effects considered here, up to terms of 
order e^(md — mu). 

We comment on this result. First, in comparison to the pion form factor 
discussed above, there is an additional term present, proportional to — m^. 
This is in agreement with the statement made before: the phase is not fixed by 
the TTTT amplitude alone, which does not contain any — rriu terms at one-loop 
order. Second, the difference to the isospin symmetric phase now grows with s 
even at one-loop order, 

^o-5o = ^^ + 0(1), s/M^^»l. (55) 

According to the discussion in section [3l one has to subtract the difference iIjq — 6q 
from the measured phase shift before a comparison with the chiral prediction can 
be performed. In figure M we display this difference in the relevant decay region, 
for i? = 37 ± 5. The width of the band reflects the uncertainty in R. It is seen 
that the isospin correction is quite substantial - well above the systematic and 
statistical uncertainties quoted for the measured phase shift [7]. 

In addition, similarly to the case of the scalar form factor, the phase-removed 
form factor for the Ke4, decays is not holomorphic at the charged threshold and, 
consequently, in the expansion of this form factor even and odd powers of the 
variable q are present, in contrast to the isospin symmetric case, see Eq. (fTOj) . If 
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Figure 9: The isospin breaking correction that must be subtracted from the phase 
shift 5 measured in K^.^ decays. The width of the band reflects the uncertainty 
in the ratio i? = 37 ± 5. 



the value of the form factor at g = is factored out, the coefficient of the term 
linear in q is the same as in the case of the scalar form factor at leading order in 
ChPT. On the other hand, this is a two-loop effect and therefore suppressed in 
magnitude. To obtain a reliable estimate of its size is beyond the scope of this 
work. 



7 TTTT scattering lengths from K^a decays 

The first large-statistics experiment to measure the phase ipQ — ipi in K^.^ decays 
has been performed by the Geneva-Saclay collaboration about thirty years ago. 
In this experiment about 30'000 K^.^ decays have been collected and analysed [4]. 
More recently the E865 experiment at Brookhaven [4] and the NA48/2 experiment 
at CERN [5] have each collected more than ten times the statistics and have made 
possible a precise extraction of the scattering length aj]. So precise in fact, that 
a proper treatment of the isospin breaking corrections as discussed in this paper 
becomes essential. 

In this section we discuss how the isospin breaking corrections infiuence the 
extraction of the scattering length, and will do this for all three data sets. We 
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Figure 10: S-wave scattering lengths plane. Contours shown correspond to one 
sigma. The ellipses do not take into account any theoretical uncertainty. 



perform two kinds of analyses: we will either leave both S-wave scattering lengths 
free and fit them to the data, or use a low-energy theorem which relates both of 
them to the scalar radius of the pion and end up with a one-parameter fit [30] . 

For the two-parameter fits we have used the parametrization of the tttt phase 
shifts corresponding to solutions of the Roy equations as functions of the two 
S-wave scattering lengths which was provided in [31]. To evaluate numerically 
the isospin breaking correction fl52l) we have used i? = 37 ± 5 and F = 86.2 ± 0.5 
MeV [32] - moreover we use the fact that tpi = 6i at one-loop order. The results 
of these fits are shown in Fig. [10] as one-sigma contours (68% probability, i.e. 

= Xmin + 2-31) for the two most recent data sets, with and without taking into 
account isospin-breaking corrections. The ellipses corresponding to the Geneva- 
Saclay data are too large in comparison and for this reason we have not shown 
them in the same figure. The figure shows that while the ellipse corresponding 
to the uncorrected NA48 data does not overlap with the theoretical prediction, 
after applying the isospin breaking correction, the two overlap completely. For the 
ellipses corresponding to the E865 data the situation is reversed: before applying 
isospin breaking corrections there is a perfect overlap, whereas after applying 
them the experimental one-sigma contour and the theory prediction barely touch. 
The figure also clearly shows that the two sets of very high statistics data are 
not in very good mutual agreement: the two one-sigma ellipses (with or without 
isospin breaking corrections) do not have any significant overlap. It has been 
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recently pointed out by Brigitte Bloch-Devaux that the origin of this tension hes 
in the data point in the last bin of the E865 data set [7], for which the evaluation 
of the barycenter may need to be revised [33]: we have verified that indeed it 
is enough to remove this point from the fit to obtain an almost perfect overlap 
between the two ellipses. 

We take into account the constraint of the pion scalar radius as follows: we 
use the numerical estimate (r^)j, = 0.61 ±0.04 fm^ [2], which implies the relation 



a, 







/(Aa°) ±0.0008, Aa° = a°~0.22, 



f(x) = -0.0444 ± 0.236 a; -0.61x^-9.9 3;% (56) 

where the error accounts for the various sources of uncertainty in the input used 
in the Roy equation solutions. In our fits we minimize the following x^'- 

_ / ag-/(Aa°) y 

V 0.0008 J ■ ^ ' 

Since the K^^ data are very little sensitive to a^, the minimum of the lies 
then always on the /(Aoq) line, but in this way we also take into account the 
uncertainty in the scalar radius relation. The results of the fits are: 

without applying isospin breaking corrections 

0.243 ± 0.037 = 2.2 Geneva-Saclay [ 5 data pts.] 

0.218 ± 0.013 x^ = 5.8 E865 [ 6 data pts.] (58) 

0.245 ± 0.007 = 9.6 NA48 [10 data pts.] , 

and after applying isospin breaking corrections 



0.222 ± 0.040 = 2.1 Geneva-Saclay 

0.195 ±0.013 x^ = 6.6 E865 (59) 

0.223 ±0.007 x^ = 11.5 NA48 . 



Averaging the latter three independent determinations yields Og = 0.217 ± 0.008 
where the error has been inflated by an S* = 1.3 factor according to the PDG 
prescription (the x^ of the average of the three fits is equal to 6.3). Repeating 
the same procedure after having removed the last data point in the E865 set 
yields Oq = 0.221 ±0.007 with x^ — 0.7, so confirming the observation of Brigitte 
Bloch-Devaux [7]. We also briefly comment on the vertex corrections at order 
aQED mentioned in section [31 They lower the value of by ~ 0.013. For reasons 
already mentioned, we do not take these corrections into account in the present 
work. 

The extraction of the scattering length Oq from the phase shift '0o ~ '0i mea- 
sured in data requires theory input in two steps: the first is the subject of 
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this article, the evaluation of the isospin breaking corrections which lead from 
i^o — i^i to 60 — 61, and the second is the step from the phase 6q — 61 to the 
scattering length. Both steps can be made quite accurately, but both are subject 
to a theory uncertainty which, though small, has to be properly accounted for 
at the level of precision currently reached by experiments. As far as step one 
is concerned, the evaluation of the isospin breaking correction given in Eq. fl52l) 
relies on an estimate of F and R as discussed above. A one-sigma change in F 
has no visible effect in the fitted value of Oq, whereas a change of R by five units 
shifts by 0.001. Our one-loop estimate of these effects is subject to higher- 
order corrections. We have evaluated these in the case of the scalar form factor 
of the pion discussed in previous sections. If we use the latter as our estimate 
of higher-order isospin breaking effects and repeat the fits, we get a value of Oq 
shifted by 0.004. An alternative quick way to estimate two-loop effects is to 
substitute F with the physical - doing this the fitted value of the scattering 
length changes by 0.003. We conclude that 0.004 is a reasonable estimate of the 
uncertainty due to higher-order isospin breaking corrections. As for step two, the 
main source of uncertainty in the calculation of how the Roy equation solution 
depend on Oq is represented by the value of the 6q phase at 0.8 GeV [31], which 
had been estimated to be (82.3 ± 3.4)°. Changing the latter by one sigma and 
repeating the fits leads to a shift of Cq by 0.004. If we sum in squares these four 
different sources of uncertainty (namely: value of F, value of R, higher-order 
isospin breaking corrections and value of 5o(0.8GeV) we end up with 

a° = 0.217 ±0.008cxp± 0.006th , (60) 

which represents the current best experimental determination of the isospin zero 
S-wave scattering length coming from KeA decays. NA48 is currently analysing 
more statistics [7], and we expect that when the updated analysis will become 
available, the experimental error will reach or even become smaller than the 
theoretical one. Other competitive determinations coming from the cusp in — > 
Sir decays [10] and from pionium lifetime [9] are also expected to be improved in 
the near future. 

8 Comparison with other work 

The isospin-breaking corrections to the Kg^ decays have been evaluated in the 
past by using various theoretical settings. We briefly review some recent articles 
where the problem has been addressed. 

In Ref. [16], a full set of isospin-breaking corrections to the Ke4 decays of 
charged kaons has been calculated at one loop in ChPT. Our result for the elas- 
tic phase shift given in Eq. (152!) agrees with the result of Refs. [16] in the absence 
of virtual photon contributions. As already mentioned in section [3l these cal- 
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culations might provide a basis for an algorithm which removes isospin-breaking 
corrections from the raw experimental data. 

In Refs. [13,14], isospin-breaking corrections to Ke4 decays have been treated 
in an approach which may be considered a merger between i^-matrix theory and 
a conventional quantum-mechanical framework used to study Coulomb interac- 
tions in the final state. Although some of the expressions in Refs. [13,14] closely 
resemble the pertinent formulae of the present work, the framework used is in- 
complete and, as far as we can see, incorrect. To mention two examples, we note 
that a relation between the amplitude considered there and the K^^ decay ma- 
trix element is not provided, and the effect caused by the quark mass difference 
T^d — f^u is not discussed at all. Moreover, according to the explicit formulae 
given in Refs. [14], the effect of isospin breaking - in the absence of Coulomb 
interactions - can be expressed through the S-wave vrvr scattering phases (in the 
isospin limit) and the pion mass difference. As we show in the present article, this 
is not correct. The authors of Refs. [13, 14] also provide a set of electromagnetic 
corrections which include both, vertex corrections and rescattering of the virtual 
pions through the exchange of Coulomb photons (the latter is a two-loop effect 
in our terminology). As already mentioned, it is not clear which part of these 
corrections was already taken into account in present analyses. We conclude that 
one should not use the results of Refs. [13, 14] in data analyses of K^^ decays. 

9 Summary and conclusions 

Data on Ke4 decays allow one to measure the difference of the S- and P-wave 
phases ipQ—ipi of a particular form factor in the matrix element of the strangeness- 
changing axial current. If isospin symmetry were exact, ipo — ipi would coincide 
with the difference of the S- and P-wave vrvr scattering phase shifts 5o ~ (Wat- 
son's theorem). We have shown that the situation changes if isospin symmetry 
is broken. In particular, 

i) The phases ipQ and 5q are not equal. Moreover, ipQ is not determined by the 
TTTT scattering amplitudes alone (e.g., the former contains the LECs which 
are absent in the latter). 

ii) The phase ipQ does not vanish at the threshold s = 4M^. 

iii) At threshold, the expansion in Eq. ( ITOi) . which is valid in the isospin sym- 
metry limit, becomes e~*'^°/o(s, s^) = cq + ciq + C2q^ + O(g^) , where the 
coefficients Cj depend on the variable s^, and similarly for the P-wave form 
factor /i (although at higher order in the chiral expansion). 

Finally, in order to determine the isospin symmetric phases from data, one 
uses Eq. f|T7|) and calculates the difference ipk — ^k, which must then be subtracted 
from the measured phases ipk- We have performed this calculation at one loop 
in ChPT for K^^ decays and have shown that the correction is substantial, well 
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above the statistical and systematic errors in present NA48/2 data. Once this is 
taken into account, the determination of the isospin zero S-wave scattering length 
from Ke4 data yields 

a[j = 0.217 ±0.008exp ± 0.006th , 
in excellent agreement with the chiral prediction. 
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A Notation 



In the main text, we use the following notation for cut complex planes: 

CL/j(a) = {z\z ^<C,z ^ (-00, 0], z ^ [a, 00)} , 
C^(a) = G C,^ ^ [a,oo)} . (A.l) 

- The region Cii^(a) denotes the complex plane, cut along the negative real 
axis, as well as along the positive real axis for 2 > a. 

- The region CR(a) denotes the complex plane, cut along the positive real axis 
for z> a. 



B Threshold behaviour of phases and form fac- 
tors 

We prove some of the statements made in section H] and start with the isospin 
symmetric case considered in subsection 14.11 The issue was discussed in the 
literature long ago [24] - we keep the presentation therefore short. 

B.l Isospin symmetric case 

The partial waves are denoted by t^. We assume that these i) are holomorphic 
in the cut plane C/ shown in Fig. IB.lt h) are real on the real axis < s < 4M^ , 
and iii) converge to continuous boundary values tfe(s) as 2 — > s. The form factors 
are boundary values of functions F^^i^z^ which are assumed to be i) holomorphic 
in the complex 2;-plane, cut at 2; = s for s > 4M^ , ii) real on the real axis for 
s < 4M^ , iii) nonzero for finite \z\. Furthermore, we assume that F^^z) converge 
to continuous boundary values as ^ — > s. 



We introduce the two-particle irreducible amplitude 

t(z) = ^4^^^. (B.2) 

It is meromorphic in C/. The poles are due to possible zeros in the denominator. 
In the elastic region, the discontinuity of t vanishes. lYiom. the Edge-of-the- 
Wedge theorem (EWT) it then follows that t is meromorphic in the cut plane 
Cii displayed in Fig. IB.li Solving for to gives 



tQ{z) = F{z)+iaG{z 

F = ^ 

1 + aH 



G = tF. (B.3) 
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Figure B.l: Regions of analyticity in the complex 2;-plane, isospin symmetric 
case. We assume that tk are holomorphic in the cut plane C/. The two-particle 
irreducible amplitude t is then meromorphic in the cut plane C//, and the phase- 
removed form factor F^. is holomorphic in the shaded region Cm. The symbol 
Zinei denotes the inelastic threshold z = 16M^. 



The quantities F, G are meromorphic in Cjj. From this representation, it follows 
that the analytic continuation of tg to the second Riemann sheet is [24] 

= ^° , (B.4) 

which shows that the zeros of the S'-matrix on the first Riemann sheet correspond 
to poles of the partial waves on the second Riemann sheet, at the same value of 
z. 

Now consider the form factor F^, and define the ratio 

m = ^. (B.5) 

which is meromorphic in Ci. Unitarity and EWT show that the elastic cut is 
absent. Therefore, one has the decomposition 

F^{z) = A{z)+iaB{z) , (B.6) 

with y4, B meromorphic in C//. Furthermore, the form factor generates a pole on 
the second Riemann sheet as well, at the same place where to does. Finally, the 
phase defined in Eq. f l26|) coincides with the isospin zero S-wave txti phase shift 
in the elastic interval, because R is real there. The phase-removed form factor 
in Eq. ^I7\ is holomorphic in the shaded region Cm displayed in Fig. fIB.ip . 
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B.2 Isospin broken case 



The threshold structure of the form factor displayed in Eq. ( !29l) can be obtained 
in a manner quite analogously to the previous subsection. One first works out 
the threshold behaviour of the scattering amplitudes. We assume that the hk i) 
are holomorphic in the region Dj displayed in Fig. IB. 21 ii) are real on the real 
axis in the interval < s < 4M^o, and ii) converge to continuous boundary values 
hk{s) as z s. The form factors are assumed to satisfy the same conditions 
as above, with M^o- 



4M^o ml Zinel. 



cut plane Dj 



inel. 



cut plane Dn 



4M^% Zinel. 



region -D/// 



Figure B.2: Regions of analyticity in the complex z-plane, isopin broken case. We 
assume that are holomorphic in the cut plane Dj. The two-particle irreducible 
amplitude T is meromorphic in the cut plane Du, and the phase-removed form 
factor Fc is holomorphic in the shaded region Dm, with the cut indicated. The 
symbol Zinei denotes the inelastic threshold z = IGM^q. 

We introduce the matrix 

compare the comments after Eq. (122]) . Consider the two-particle irreducible 
amplitudes 

f={l + iTpy^T, (B.8) 

where the matrix T is defined in Eq. ( 12T|) . Unitarity and EWT show that the cut 
4M^o < s < 16M^o is absent for T, from where we conclude that the matrix T 
has the threshold structure 

T = A' + iaB' + iaoC + aa^D' (B.9) 
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with matrices A', B', C, D' whose entries are meromorphic in Du. Next, one 
forms 



Rc 

Rq 



(B.IO) 



where F is the matrix in Eq. ( 12T|) . and uses again unitarity and EWT to show 
that -Rc,o are meromorphic in Djj. Solving for F gives 



Fk = Ak + iaBk + iaoCk + cra^Dk ; k = c,0. 



(B.ll) 



with Ak,-- - ,Dk meromorphic in Djj. Evaluating the phases 6k according to 
Eq. (130|) shows that these depend on the ratio Rq/Rc- On the other hand, this 
ratio depends on Fc/Fq, which is unity in the isospin limit, and different otherwise. 
In subsection 15.21 we show that, in case of ChPT, Fc/Fq contains LECs that do 
not occur in the tttt scattering amplitude. Hence the phases 6k are not determined 
by the tttt amplitude alone in the case of broken isospin. From fl30l) . it is seen that 
the phase 6^ does not vanish at the threshold s = 4:M^, because ao(4M^) 7^ 0. 
Finally, from Eq. fIB.lip . one finds for the phase-removed form factor the result 



Fc = {Ac + aaoDc 



Ac + aa^D^ 



(B.12) 



It is holomorphic in the cut region Dm displayed in Fig. IB.2[ The cut is due to 
the factor aao. 



C Non-relativistic effective Lagrangian frame- 
work 

The non-relativistic Lagrangian C = Ct,t, + Co, which is used for calculation of 
the scalar form factor, consists of a part describing quartic pion-pion interactions, 
and a part corresponding to the interaction of the pion pair with the external 
current 0{x). The framework is described in detail in our recent works [25-27] 
and will not be repeated here. For example, the Lagrangian Cj^tt is given in Eq. 
(4) of Ref. [25]. The polynomials di{z), introduced in section [5AI are given by 

327r4 = C+„ + (z - 4M^)D+_ , 

327r4 = C, + {z-4M^)D,, 

327rdo = Coo + {z-AM^o)Doo. (C.13) 

Here we use the same notations for the 4-pion couplings C, D as in Refs. [25-27]. 
The matching of these couplings to the effective-range expansion parameters in 
the vrvr scattering amplitudes is described in these references as well. 
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The part of the Lagrangian which is responsible for the interaction with the 
external source takes the form 

Co = O (^-fi<^l<^l + f2{W±^lW±<^l + \/<^lV<^l- M^<!?1<^1) 

+ y ^Wo - f (Wo^lWo^l + V$SV$S - Ae,$;$;)) + h.c. (C.14) 

where $-|-,$o denote the non-relativistic charged and neutral pion fields, W± = 
\/ Ml — A and VFq = y^M^o ~ A are the pertinent differential operators, and 
/i; ■ ■ ■ /4 stand for the various non-relativistic couplings. The polynomials /c, /o 
which are used in the main text are defined as 

fc = h-\{z- ^Ml)h , h = h-\{z- ^Ml,)f, . (C.15) 

In order to calculate the scalar form factors within the non-relativistic framework, 
one evaluates the transition amplitude between the non-relativistic two-pion 
states and the vacuum in the presence of the external source The form 

factors are defined by 

(7r+(pi)7r-(p2);out|0;in)^^ = i j d^xe'^P'+P'^''0{x) Fc{s) + 0{0'^) , 

(7r°(pi)7r°(p2);out|0;in)^^ = t J d^xe'^P'+P'^^'Oix) Fo{s) + 0{0^) , 

(C.16) 

where s = (pi +^2)^- 

The behaviour of both the relativistic and non-relativistic form factors at 
small momenta is the same and is displayed in Eq. fl29|) . The polynomials fk 
can be determined by performing a matching of the regular part of the form 
factors, evaluated in both theories. 
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